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Pure-gauge SU(2) on large lattices 

UKQCD collaboration, presented by P.W. Stephenson** 

'^Department of Applied Mathematics and Theoretical Physics, Chadwick Tower, 
University of Liverpool, P.O. Box 147, Liverpool L69 3BX, United Kingdom 
Email: pws@s-a.amtp.liv.ac.uk 

We have calculated the qq potential over a wide range of lattice separations on a 48^ ■ 56 lattice at f3 — 2.85. We 
are able to investigate both long-range and (by correcting for the effects of lack of rotational invariance) short- 
range potentials. From the former we estimate the string tension and from the latter we are able to investigate 
asymptotic scaling. Fitting to this enables us to give an estimate for the A-parameter of pure-gauge SU(2). 



1. Introduction 

We discuss recent results for the potential be- 
tween (hypothetical) static quarks calculated on 
the 64-node Meiko Computing Surface at Edin- 
burgh University. Our lattice size, 48^ • 56, is to 
our knowledge the largest so far used in a lattice 
gauge theory simulation. Thus we are in a good 
position to examine the approach to the contin- 
uum limit and comment on the remaining effects 
of finite lattice spacing a. A fuller write-up is 
available as a preprint Q, so after briefly men- 
tioning our data and techniques we proceed di- 
rectly to a discussion of the results. 



2. Monte Carlo data 

We have chosen (3 = 2.85 based on previous ex- 
perience to make our lattice physically as large as 
necessary. It is known that the potential is less 
sensitive to finite-volume effects than, for exam- 
ple, the glueball spectrum but we do not wish to 
rely on this. We also have results for glueballs 
which will not be discussed here. We have per- 
formed a total of over lO^'^ link updates on this 
lattice. 

Our technique for extracting the potential Vqq 
is standard. The effective potential in lattice 
units is extracted from R x T Wilson loops as 
Veff{R, T)a = ln(l^(i?, T)/W{R, T- 1)); as T be- 
comes large, Ves{R,T) tends to the ground state 
potential Vqq{R) which is what we need. To im- 
prove our overlap onto this state we have used two 



popular 'fuzzing' schemes (|,^; the 'wave func- 
tions' generated have been used as the basis for 
a variational analysis. The resulting overlap ap- 
proaches unity in many cases. 

We are faced as always with the problem of 
deciding how to extract Vqq{R) from T4fi(i?,T). 
We find that the results at T = 5 agree within 
statistical errors with those at T = 4 and use 
these as our main results. We have also performed 
an extrapolation to T = oo by assuming that the 
the signal is contaminated only by one excited 
state. This assumption gives a lower bound on 
the potential. In this case we are less in control 
of statistical errors; furthermore, we expect this 
procedure to be less effective for larger R where 
the first excited state is closer to the ground state. 
Thus we use the extrapolated potentials mainly 
as a way of gauging our systematic errors. 



3. Fits to data 

We have performed two separate fits for large 
and for small separations R. In fact, we have 
used different data analysis techniques in the 
two cases: the large-separation analysis required 
larger statistics but we reduced the total amount 
of data produced by blocking the lattice by a fac- 
tor of two in each spatial direction. In the other 
case we needed the full lattice volume but lower 
statistics for accurate results. 
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3.1. Large R 

Here our major object is to investigate the 
string tension K. We fit to the form 

ViR)=C-^+KR (1) 

for R > 6a. In fact, the Coufomb term is modified 
to include corrections (from first order perturba- 
tion theory) for lack of rotational invariance, but 
the effect in this case is negligible. We are not in- 
terested in the constant C since the ground state 
energy is in large measure an artefact of the cut- 
off. Our results for both the T = 5 and extrapo- 
lated potentials are shown in table |l 



6G{R) is the one-gluon correction to the Coulomb 
term, so that if / were unity the two terms in- 
volving A would combine to give the first-order- 
corrected Coulomb term which was used in the 
previous fit. By allowing / to vary we account 
for the following observation: that the dominant 
effect of higher-order corrections is to rescale the 
first-order term. 



Table 2 

Result of fits for R > a (with K fixed) 
A f B/a xVdx^ 

0.238(7) 0.63(3) 0.050(6) 9.7/14" 



Table 1 



Result of fits for R> 6a 





E 


Ka^ 


X'^/d.o.f. 


V{R,5) 


0.247(7) 


0.00401(8) 


15.2/7 


Voo{R) 


0.247(14) 


0.00363(17) 


8.2/7 



As already mentioned, we use T4o to estab- 
lish a systematic error. Hence the value we use 
for the string tension is from y(i?, 5): Ka^ = 
0.00401(8)(38). 

3.2. Small R 

Here our fit includes all separations R > a. 
Furthermore, we also include off-axis potentials 
— that is, those with separation vectors such as 
(1, 1,0), (1, 1, 1), and so forth. This makes some 
means of controlling lattice artefacts imperative. 
Our fit in this case has the form 

V{R)^C-^ + ^ + KR-Af5G{R), (2) 

in which K is fixed from the preceding fit; A, B 
and / are the free parameters. 

There are two new elements. In this range of 
separations we are hoping to be able to see the 
coupling run; the term in i3 is a simple-minded 
way of correcting the Coulomb term to account 
for this over a small range of R. 

Both A/R and BjR? terms describe effects 
which are smooth and monotonic in R. By con- 
trast, the final term describes lattice artefacts: 



The of the resulting fit (table H) justifies our 
use of this ansatz to correct for lack of rotational 
invariance. However, it is difficult to estimate the 
systematic error one should attach to this. We 
have arbitrarily assigned 10% of the correction 
which has been applied. 

4. Scaling 

The running coupling is conventionally de- 
fined ^] in terms of Vqq as 

( Ri+R2 \ 4 ^ ^ V,{R{) - K(i^2) 
a,,\-^)^^-R,R, , (3) 

in which we use the potentials Vc(-R) with the 
lattice-artefact correction applied. This is a con- 
venient quantity to use in investigating scaling as 
it is dimensionless. 

We therefore compare these results with those 
obtained previously Q on 32"^ lattices at /3 = 2.7, 
plotting oiqq against R and rescaling only the 
horizontal axis. That the agreement is entirely 
comfortable, with no indication of scaling viola- 
tion, will become clear when we discuss figure |l| 
which imposes the even stronger condition that 
the rescaling factor is \fK. 

For now, we show in table ^ various val- 
ues for the ratios of the lattice spacings 
a(/3=2.7)/a(/3=2.85). The first value is obtained 
by a full fit to the Monte Carlo results, including 
an estimate of the systematic error. The last is 
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Figure 1. The running couphng versus a length scale 



Table 3 



The ratio a(2.7)/a(2.85) in various calculations 



UKQCD 


scheme 


Liischer 
et al. @ 


Asymptotic 
(two loops) 


1.60(6) 


1.53 


1.48 


1.46 



the naive expectation of 2-loop perturbation the- 
ory using the bare couplings. It is clear there is 
some discrepancy. The other two values require 
some comment. 

It is useful to use instead of the bare cou- 
pling something more closely related to physi- 
cal processes. One popular choice is the /3e- 
scheme , where the plaquette expectation value 
(C/pi) is used as the basis for calculating the cou- 
pling. Here we use the simplest variant where 
Pe = C2/(l — {Up\)), in which C2 is a perturba- 
tive coefficient. Using even this simple correction 
changes the a-ratio by about half of the discrep- 
ancy between Monte Carlo and nai've perturba- 
tion theory, giving some confidence our results 



are in the right region. 

The remaining estimate is taken from the fit 
in figure 3 of ref. Q. In this paper (see the 
talk by Martin Liischer at this conference) the 
authors calculate the running coupling by a re- 
cursive finite-size-scaling technique. The value 
shown is in terms of the bare coupling with lattice 
sizes determined at fixed physical coupling. The 
result quoted, close to the two-loop value, is not 
central to ref. ||] ; it is indicative of the fact that 
the starting point for that paper is in small (and 
therefore perturbative) volumes. 

We have also compared the new results with 
those on 24"* lattices at /? = 2.4. Even here we 
find excellent agreement, showing scaling over a 
range of lattice spacings a(/3=2.4)/a(/3=2.85) = 
4.12(2), which is certainly impressive. 

5. Running coupling 

Equation ^ defines a running coupling in the 
Aij-scheme; A^j — 1.048A-pg so these schemes 
are perturbatively close. We plot Uqq against the 
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dimensionless Ry K which gives a physical length 
scale. This plot is shown in figure |l|. 

In the figure we show results for both j3 = 2.85 
(O) and (3 = 2.7 (□). For clarity only the larger 
of the two error bars (statistical and systematic) 
is shown here. (Both sets appear in ref. 0.) The 
effect of running is clear and the curve appears to 
be universal. 

The solid lines arc plots of aqq from perturba- 
tion theory: 
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(4) 



4^ (6o ln(Afl^i?)-2 + (61/60) In ln(Afl.i?)-2) 

with the usual perturbative coefficients. The up- 
per line corresponds to ah-n = 0.044 and the 
lower to ah-n = 0.038; note that our results in- 
clude increasing contributions from the linear po- 
tential as R increases which will tend to make 
our results larger than the perturbative estimate. 
From these limits our estimate of the scale pa- 
rameter is aAij = 0.041(3). In terms of the 
string tension this is V^/A^ = 1.54(15) or 
\/K/Ajjg — 1.62(14); using the (in this case 
entirely unjustifiable) value = 0.44 GeV we 
get A;^ = 272(24) MeV. (In these units our 
shortest-distance point, at i? = V^a, would corre- 
spond to a momentum scale q = \/ Rk, 4.4 GeV.) 

The remaining points with horizontal error bars 
are converted from ref. We have used the 
value of the string tension calculated here to set 
the scale, however the error in the string tension 
is not included. (This is consistent with the other 
points in the figure.) The region of interest is at 
the lower end; the error bars just fail to overlap 
our lines for A^, so the conclusions are not en- 
tirely clear. If our estimate of systematic errors 
is to be believed there appears to be some dis- 
crepancy. 

6. Concluding remarks 

It has recently become clear [Q that attempts 
to compare Monte Carlo results from the lattice 
with perturbation theory using the bare coupling 
are doomed. Here we have presented confirma- 
tory evidence which further suggests that when 
one uses more physical quantities, the usual Wil- 



son theory for SU(2) without fermions is able to 
produce results entirely in accord with those in 
the continuum. Our bare coupling is not able to 
give us agreement with perturbation theory, but 
provided we fit to the A-parameter ourselves we 
see a running coupling entirely in agreement with 
the asymptotic two-loop result. In addition we 
see scaling of Vqq results between lattices differ- 
ing by a factor of four in lattice spacing. 

In ref. |^ these new results are displayed in 
the I3e scheme referred to above and it is clearly 
seen that the approach to asymptotic behaviour 
is much smoother when the bare coupling is not 
used, enabling an extrapolated value of the ra- 
tio ^/K/Ajjg to be extracted. Their value is 
1.79(12), slightly higher than our unextrapolated 
number which is 1.62(14). 

Overall, we believe that pure-gauge lattice 
SU(2) is already able to give results which ac- 
curately reflect continuum dynamics. 



S.P. Booth et. 
large lattices. 
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